Optimal measurement precision of a nonlinear interferometer 
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We study the best attainable measurement precision when a double-well trap with bosons inside 
acts as an interferometer to measure the energy difference of the atoms on the two sides of the trap. 
We introduce time independent perturbation theory as the main tool in both analytical arguments 
and numerical computations. Nonlinearity from atom-atom interactions will not indirectly allow the 
interferometer to beat the Heisenberg limit, but in many regimes of the operation the Heisenberg 
limit scaling of measurement precision is preserved in spite of added tunneling of the atoms and 
atom-atom interactions, often even with the optimal prefactor. 
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I. INTRODUCTION 

Paraphrased for optical or atomic interferometry, the 
Heisenberg limit [l| states that the best possible achiev- 
able uncertainty of phase measurements is inversely pro- 
portional to the number of bosons N. The Heisen- 
berg limit, with continuous refinement and sharpening 
of the concept 0, Q , is one of the enduring paradigms 
in quantum metrology. Nevertheless, there is nothing sa- 
cred about 1/N scaling; in nonlinear schemes where the 
quantity to be measured couples to a fc-body operator 
the measurement uncertainty could scale as 1/N k 0, @. 
This is well and fine if the aim is the best possible mea- 
surement of, say, the atom-atom scattering length [fl-[H , 
which in fact enters the time evolution in conjunction 
with a two-body operator. However, such a 1/N 2 scal- 
ing docs not directly help if you insist on a measure- 
ment of a quantity that does couple to the first power of 
boson number. The nonlinear (in light intensity) Fara- 
day rotation, a possible method to measure a magnetic 
field, demonstrably allows for improved noise properties 
that originate from the nonlinearity [H , but virtually all 
of practical interferometry is based on one-particle cou- 
pling. 

The impetus to the present work is the question of 
what happens to the Heisenberg limit, prefactor or scal- 
ing, when the time evolution of a probe depends on the 
quantity to be measured and also has other components 
that modify the evolution, particularly nonlinearity. Our 
analysis takes place in two steps. 

In Sec. |H]we address the general case. Schematically, 
the dynamics that converts the quantity to be measured 
9 into a state change of the probe has two parts, one 
described by a Hamiltonian OG, where G could be called 
generator of 9 translations, and an additive part K that 
does not depend on 9. We show how time independent 
perturbation theory may be harnessed to analyze the best 
possible precision of the measurements of the quantity 9, 
and subsequently argue that the addition of the term K, 
any term K, cannot improve the precision @, @. 

Next we formulate a tangible example, measurement 
of the energy difference of an atom in two potential wells 
in the presence of both tunneling of the atoms between 



the wells and atom-atom interactions. The same physi- 
cal system was studied for similar aims in Ref. 0, al- 
though from a more practice oriented standpoint. The 
optical analog would be a device that allows tunneling 
(exchange) of photons between the two arms of the in- 
terferometer and also has nonlinear phase shifts, as per 
the Kerr effect. 

We state the example in Sec. IIII1 and analyze the 
achievable limit of precisions in Sec. IIVI as a function of 
the parameters of the system. Some results, such as how 
the interplay of tunneling and atom-atom interactions af- 
fect the precision, should be educational per se, but more 
to our point of principle, we have here explicit examples 
in which the pieces added to the time evolution of the 
probe reduce the best possible attainable measurement 
precision. For the most part the Heisenberg limit 1/iV 
scaling is still retained, but with a reduced prefactor. A 
few remarks in Sec. |V] wrap up the paper. 



II. GENERIC MEASUREMENT SCHEME 

We start with the general scheme of quantum measure- 
ments: Prepare an initial state for the "probe" , have the 
probe evolve under some quantum mechanical law that 
depends on the quantity to be measured, and finally in- 
fer the value of the quantity from measurements on the 
probe. The best possible precision in the setup we are 
considering is achieved with a pure initial state 
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we take the measurement to start with a pure state 
Moreover, we assume that the quantum mechanical evo- 
lution is generated by a hermitian Hamiltonian. A pure 
state then remains pure during the evolution. We have 
a unitary mapping of the initial state |^o) to the state 
\tp{9)} that depends on the parameter 9 to be measured, 
and we may write 

^ K{8) IVo) , 



\i>) = \m) 

where K(9) is Hermitian. Let us also define 



W) 



W(6)) = ^\W)) 



(1) 



(2) 



The criterion we use for the precision of the measure- 
ment is the variance, the square of the standard devia- 
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tion. A wcll-dcfincd procedure exists to find the smallest 
value of the measurement uncertainty over all possible 
measurements, given the initial state and the actual value 
of the parameter 9 (!), and it even describes a measure- 
ment that could be carried out (in principle) to reach 
the minimum uncertainty [l^ - fl^ ]. The variance is ex- 
pressed as the inverse of quantum Fisher information F, 
which in turn is found from the symmetric logarithmic 
derivative of the density operator A. A straightforward 
calculation [Tij using the norm conservation 



<V#'> + <v'M = o 



(3) 



shows that the symmetric logarithmic derivate, the 
Fisher information, and the standard deviation of the 
measurement results equal 



A = 2(|/) M + |V) (V'l), 

F = Tr(pA 2 )=4[(^ / )-KV'IV'')| S 



ere 



1 

7f 



(5) 
(6) 



It pays to notice that everything we say in the present 
Sec. |IT] is completely general. Every measurement that 
starts from a pure state and is based on Hamiltonian 
time evolution where the quantity to be measured is a 
parameter in the Hamiltonian is covered. 



A. Perturbation theory 

The technical innovation here is to employ time inde- 
pendent perturbation theory in both numerical and ana- 
lytical work. We explain presently 

The equality 



IV-') 



do 



-iK{8) 



(7) 



guides us to examine the derivative of the evolution op- 
erator e~ lK ^ with respect to the parameter 8. For an 
infinitesimally small dd we have 



d0 — e - lK ( e ) ~ e -iK(8+d8) _ e -iK(8) 

d9 

~ p -i(K+d6K') _ -IK 



(8) 



where we have dropped a few explicit arguments 6 and 
defined 



K 1 = K'{6) = -K{9) 



(9) 



Suppose we know the eigenvalues A n and eigenvectors 
|n) of K(9), so that we have the spectral representation 



(10) 



One way to proceed with Eq. © is to try and find the 
similar spectral representation for the operator K + d6 K' 



for an asymptotically small d9, eigenvalues and eigenvec- 
tors of the operator K perturbed by the "small" operator 
dO K 1 . This is evidently an exercise in the usual time in- 
dependent perturbation theory. 

Let us first assume that the eigenvalues A n of the op- 
erator K arc nondegencratc. To the leading nontrivial 
order in d9 the eigenvalues and eigenvectors of K + d8 K' 
are then of the form 



A n =A n + d0£„, \tp n ) = \n) + dO \(f> n ) , (11) 



with 



C n = (n\K'\n), \ef> n )=^\m) H K'\n) ^ 



rn^n 



A n — A, ; 



If there are degeneracies, we are dealing with a combina- 
tion of degenerate and nondegenerate perturbation the- 
ory. We then choose the eigenstates of K so that in any 
degenerate manifold they are also eigenstates of K', and 
simply drop the seemingly divergent terms in the sum 
in (|12[) . To the leading nontrivial order in dO we then 
have 

e - i(K+MK ') ^e-iA^^ 

n 

n 

+G^]Te- 4A " (|^ n )(n| - i£„|n)H + |n)<0„|) , (13) 

n 

which gives 
d 



de 



(14) 

The model that thoroughly permeates quantum 
metrology literature stipulates 



K{6) = 6G, 



(15) 



where the hermitian G could be called generator of 9 
translations. In this case we find 



W) 



-iG\ 



-iGe 



-i9G 



(16) 



both directly from Eq. and also indirectly from per- 
turbation theory using Eqs. ©, (HH), and The 
Fisher information is then nothing but four times the 
variance of the generator G in the state \4>o) (l2l4Iij . 

However, in all but the simplest examples like this, 
one would have to carry out the analysis numerically. 
Our method is to find the quantum Fisher information 
by combining Eqs. ©, 0, (JT4J), and (fT2]). In numer- 
ical computations perturbation theory circumvents the 
need to take any derivatives numerically. A combina- 
tion of non-degenerate and degenerate perturbation the- 
ory would make the computations tedious, so in numer- 
ical analysis we ordinarily eliminate the degeneracies by 
adding a tiny perturbation to K to break the symmetry 
that causes the degeneracies. 
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B. Limit of precision 

We add one more layer of optimization and find the 
initial state \4>o) that produces the smallest possible mea- 
surement uncertainty, i.e., maximum Fisher information. 
In other words, we maximize F(\ipo}) with respect to the 
initial state IV'o)- There are special cases when this can 
be done analytically, too. For instance, in the model (Tl5|) 
the maximal Fisher information and a corresponding ini- 
tial state are 

Fm = (9M - g m )\ - TfdffM) + \g m )) , (17) 

where <?m and g m are the largest and the smallest eigen- 
value of the generator G, and ip is an arbitrary relative 
phase between the respective eigenstates. 

As general considerations go, we first write the deriva- 
tive state from Eqs. (TTJ]), (P, and ([11]) in the form 



W) 



do 

iL\ 



-iK{8) 



dO 



e -iK{9) e iK(9) 



(18) 



where we have defined 



l 1 f A "'""' ) \ m ) ( m \ K> |n) (n\ . (19) 

-i(A„ - A m ) 



In the formally singular term with m = n the ratio is 
interpreted to have the value as appropriate for the limit 
A„ — > A m , namely 1. The hcrmitian operator L could be 
called local or instantaneous generator of 9 translations 
at the given value of 9. 

As far as it comes to the largest attainable value of the 
Fisher information, we may just as well optimize with 
respect to the state \ip), the image of the initial state 
l^o) m a norm preserving and one-to-one unitary map- 
ping e~ iK . A comparison of Eqs. ([16]), ([H]), and (fit]) 
immediately shows that the maximal Fisher information 



F, 



M 



M 



(20) 



where £m and £ m are the largest and the smallest eigen- 
value of the operator L. An optimal input state is ob- 
tained by inverting the transformation e~ lK ( e \ or 



IV'o) 



7T 



(K 



Ml 



(21) 



C. Effect of added couplings on precision 

Suppose we start with the completely solvable 
model (|15[) . and ask how the measurement precision is 
affected if additional terms need to be considered in the 
dynamics of the probe. We write 



K = 9G + K, 



(22) 



where K stands for the added dynamics that is indepen- 
dent of the parameter 9, so we have K' = G. Conversely, 
locally, around a given value of the parameter 0q, an arbi- 
trary (differentiable) K{9) can always be written in this 
form by choosing G = K'(9 ) and K = K{9 Q ) - 9 a G. 
From Eq. ([TO]) the local generator L may be written 

L = j dxY^ \m) e~ ixXm (m\ G \n) e ixX " (n| (23) 

= j dxL{x), (24) 
Jo 
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with the definitions 



L{x)=U{x)GU\x), U{x) = Y J 



(25) 

U{x) is unitary, so the spectra of all L{x) are the same 
as the spectrum of G. But now, Eq. ([24]) can be viewed 
as a linear combination of a large number of operators 
L(x) with positive coefficients that sum up to one. It is 
then easy to see from the variational principle that the 
largest eigenvalue iu of L can be at most as large as 
the largest eigenvalue among the operators L(x), i.e., the 
largest eigenvalue gu of G. Specifically, let \£m) be a nor- 
malized eigenvector belonging to the largest eigenvalue of 
L, then we have 

Im = (hi\ L \£ M ) = f dx (£ M \L(x) \£ m ) 
Jo 

r-l 

< / dx g M = g M ■ (26) 



For the smallest eigenvalue we have similarly g m < £ m . 

Therefore the maximal Fisher information (£m — £ m ) 2 
in the amended measurement governed by the local gen- 
erator L is at most equal to the maximal Fisher informa- 
tion (gM — gm) 2 for the original generator G. The added 
evolution K cannot improve the best attainable measure- 
ment precision. On the other hand, as we will see below, 
the range of the eigenvalues of L may be narrower than 
the range of the eigenvalues of G, so that the attainable 
measurement precision may decrease. 

Other authors give similar results from more [f| or 
less related arguments. Cast in terms of our present 
development, Ref. 0] in fact claims that the attainable 
measurement precision remains unchanged. The error, 
again paraphrased for the present rendition of the math- 
ematics, is in the implicit assumption that all of the op- 
erators L(x) are the same, which in general does not hold 
true even if their spectra are the same. 



III. EXPLICIT MODEL 

Our example is about a Bose-Einstein condensate in a 
double- well trap, as described by the two-mode Hamilto- 
nian (l5| or the two-site Hubbard model put in the form 



H 



-tJ x + eJ z + UJf 



(27) 
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Given the boson annihilation operators for the left and 
right halves of the potential well and a^, we have three 
operators that obey the angular momentum algebra, 



Jx = \{ a \ a R + a R a L), 

J y = m( a L a R- a R a L), 
= \i a \ a L - a R a R ) . 



(28) 
(29) 
(30) 



The parameter r is the tunneling amplitude for the 
atoms from one site to the other, e is the energy differ- 
ence for an atom in the two sites, and U represents the 
strength of the atom-atom interactions. The relevant bo- 
son states are spanned by \til, tir), where ni and tir are 
the numbers of the atoms in the left and right traps, or 
| Jm) = \J + m,J — m) that are eigenstates of the an- 
gular momentum with the components (|28|) - (|30|) for the 
total angular momentum J and its z component m. For 
a fixed number N — 2 J atoms, the state space has the 
dimension N + 1 . 



As has been known for quite a while 15|, the validity 
of the two-mode model is no longer guaranteed when the 
atom-atom interaction per particle becomes comparable 
to the energy difference between the one-particle eigen- 
states in one or the other of the potential wells. This is 
the case with experiments in the limit of the Thomas- 
Fermi approximation. However, the atom-atom interac- 
tion parameter U and the validity of the two-mode model 
may, in principle, be controlled independently by adjust- 
ing both the trapping frequency of the potential wells 
(say, by adjusting the intensity of the trap lasers) and 
the atom-atom scattering length (say, by making use of 
a Feshbach resonance). In what follows, the two-mode 
approximation is always assumed valid. 

We take it that e is the parameter to be determined. 
Time evolution according to the Hamiltonian (|2"7|) corre- 
sponds to the generator of measurement results according 
to K = Ht/H, but we absorb the factor t/H into the defi- 
nition of the parameters and make no difference between 
the operators H and K. The measurement is about a 
dimensionless evolution phase, with e standing for what 
we denoted by 9 in Sec. |TTJ 

Such a measurement is to an extent a standard task. 
In the usual optical Mach-Zehnder interferometer e would 
be proportional to the phase difference of light incurred 
between the two arms of the interferometer, and likewise 
in atom interferometer setups. If we momentarily ignore 
both the tunneling between the two sites and the inter- 
actions between the atoms, we have the usual situation 
for an atomic or optical interferometer with 



K = eJ z 



(31) 



Given N bosons, a comparison with (fT5)) and (JT7J) shows 
that the maximum Fisher information and the corre- 
sponding input state are 



With eg = 0, we manifestly have what is called the NOON 
state It gives the best possible measurement preci- 
sion with the standard deviation a e = 1/V-Fm = 1/-/V, 
the usual Heisenberg limit. 

The assignment now is to find out how the best pos- 
sible measurement precision is affected when tunneling 
and nonlinear boson-boson interactions are included. In 
fact, the added nonlincarity oc J 2 that commutes with 
the generator J z in itself has no effect on the achievable 
measurement precision. The situation becomes nontriv- 
ial only when site-to-site tunneling proportional to the 
noncommuting operator J x is also involved. 

Even though the parameters r, e and U in the Hamilto- 
nian have the dimension of energy and one of them could 
be arbitrarily picked as the unit of energy, because of the 
time evolution they get multiplied by t/H and are ren- 
dered dimensionless. This means that we already have 
three independent parameters. Moreover, on top of the 
double optimization over both the measurement process 
and the initial state as in Eqs. (0 and ([20]) we have an- 
other parameter to consider, the total number of bosons 
N. In that regard, in the special case with r = the 
maximum Fisher information is N 2 . Second, given N 
atoms, the effects of particle number in the Hamiltonian 
scale with various parameters approximately as tN , eN , 
and UN 2 . To make the physics as invariant as possible 
with respect to the number of atoms, it pays to compare 
e and NU . In short, we will express the results in terms 
of the scaled variables 



/ = F/N 2 , u = NU. 



(33) 



F 



M 



N 



\^u) = ^(\N,0) + e^\0,N)). (32) 



The Heisenberg limit is given by / = 1. 

All numerical computations were done using Mathe- 
matica. To avoid the issue of degenerate versus non- 
degenerate perturbation theory, we advise to avoid pa- 
rameter values precisely equal to zero. Even very small 
nonzero values that are in practice equal to zero eliminate 
the problems, presumably by reducing the symmetry of 
the Hamiltonian. 



IV. RESULTS 

Let us first ignore atom-atom interactions and write 

H = T(sin (f) J x + cos (f) J z ); (34) 
T = Vr 2 + e 2 , cos^ = e/T, sm<j) = -t/T. (35) 

The Hamiltonian (|34[) is proportional to the component 
of the (fictitious) angular momentum in the direction 
cos <p e z + sin <p e x leaning an angle 4> away from the z 
axis. The corresponding angular momentum eigenstates 
| Jm) j are eigenstates of the Hamiltonian, and the eigen- 
values are A,„ = m.T '. 

Computation of the local generator L from Eq. (fl9|) 
is a tedious affair involving the matrices that govern the 
transformations of the eigenstates of angular momentum 
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FIG. 1. (Color online) 3D and contour plots of the scaled 
maximal Fisher information ju = Fm/N 2 as a function of 
the tunneling amplitude r and scaled atom-atom interaction 
u — NU, given the atom number N = 2 and energy difference 
e = 1. 

under rotations (l7j , but we may produce a simple analyt- 
ical example in the limit T — > oo. Then the eigenvalues 
of the Hamiltonian A m tend to infinity, in Eq. (|19[) only 
the diagonal elements with m = n survive, and we have 
the matrix elements of the generator of local translations 

^( Jm\ L | Jm% = <5 mm / ^( Jm| J z | Jm% . (36) 

The expectation value of a component of an angular 
momentum in an eigenstatc of the component of the an- 
gular momentum tilted by an angle <f> is in fact governed 
by the classical projection, so we have the diagonal ele- 
ments of the matrix L, and at the same time its eigenval- 
ues £ m , in the form £ rn = mcos0. The maximal Fisher 
information is therefore Fm = J 2 / 4 cos 2 <j> = N 2 cos 2 </>, 
and the measurement uncertainty increases by a factor 
of 1/| cos 0| compared to the case without tunneling. In 
particular, for small values of the energy difference e the 
tunneling reduces the measurement precision by a factor 



of~|e/r|. 

We next re-instate atom-atom interactions. We begin 
in Fig. [T] with the case N = 2. Panel 1(a) is a 3D plot of 
the scaled Fisher information /m maximized over all in- 
put states as a function of the tunneling and atom-atom 
interaction parameters r and u, with e = 1. We have 
plotted only nonncgativc values of u. meaning, repulsive 
atom-atom interactions. The Fisher information is not 
exactly an even function of u, but the difference between 
u and —u for the purposes of these drawings is so small 
that it would be barely discernible in Fig. [TJ or in Fig. [2] 
below. As expected, we have the Heisenberg limit all 
along the r = axis. Increasing the tunneling amplitude 
t decreases the maximum Fisher information, while in- 
creasing atom-atom interactions counteracts the effect of 
tunneling [Toj . 

As a matter of fact, in the formal limit \u\ — > oo the 
cigenstatcs of the Hamiltonian approach the eigenstates 
of J 2 . We will eventually have doubly degenerate mani- 
folds made of states | Jm) and | J — m) , except for m = 
when there is no degeneracy. Perturbation theory re- 
quires that we diagonalize J z within each of these man- 
ifolds, whereupon we need to pick the eigenstates of the 
Hamiltonian to be eigenstates of J z itself. Therefore, in 
the case |u| — > oo Eq. (|19l) shows that the local generator 
equals the original generator, L = J z , and the Heisenberg 
limit ensues. 

Although we generally say little about the optimal in- 
put state, we mention here the special case u — > — oo. 
It is known [l8| that in this limit, and with e = 0, the 
zero-temperature ground state is a "Schrodingcr cat" or 
NOON state of the form (|3"2")l that maximizes the Fisher 
information. The practical complications are probably 
substantial, but in principle the right input state comes 
almost for free. 

Figure [5] gives similar contour plots as Fig. [TJ but for 
different atom numbers. We firstly have N = 8, 16, and 
32. The contours are all at the same values of Jm- It 
is obvious that in the limit N — > oo the scaled maxi- 
mum Fisher information = F/N 2 as a function of 
the parameters r and u = NU converges to a universal 
function /m(t, u,e,N — > oo). While the Heisenberg scal- 
ing of measurement uncertainty is distinctly a quantum 
effect, this is also the limit in which one expects that 
the semiclassical approximation becomes valid: Barring 
pathological quantum states of the system such as the 
Schrodinger cat (sic!), 6l and bp may be treated as clas- 
sical variables with certain Poisson brackets instead of 
quantum operators. 

We have investigated the convergence with increas- 
ing iV for a few fixed values of r and u. It appears 
that for large N the deviation of Jm from the N — > oo 
limit scales as TV -1 . Based on this observation, we have 
done a Richardson extrapolation [l9| to the limit func- 
tion fj\i(T 7 u,e 7 N — > oo), and plot it in Fig. [2] labeled as 
N = oo. The absolute error of the limit function is con- 
servatively estimated to be less than 10 -3 in the whole 
plotted range, which is far below the resolution of the 
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FIG. 2. (Color online) Contour plots of /m as a function of 
r and u with a fixed e = 1 for TV = 8, 16 and 32, as marked 
on the panels. Also shown is an extrapolation to the limit 
TV — > oo denoted by TV = oo. The maximum value of /m is 
attained at r = 0, and equals /m = 1. The spacing between 
the contour lines equals 0.1. 



plot. The limit function is similar to the N = 2 function 
/m (t, u, e, 2), and the preceding qualitative discussions 
of this case still apply. 

So far we have set the value of the quantity to be mea- 
sured as e = 1; next we address the variation of Fisher 
information with the parameters e itself. One conceiv- 
able reason is that the parameters e, r and u all scale 
simultaneously with the interaction time. Accordingly, 
for the time being we write e — > xe with a scaling factor 
x, and similarly for r, u. 

In the limit x — > one finds that K -> 0, e~ iK -> 1 - 
iK, \ip'} —> —iJ z \ip), and we are back to the Hciscnberg 
limit. However, the measured value of e scales to zero as 
well, and the relative accuracy becomes poor. 

To address the opposite limit x — > oo we note that K 
and its eigenvalues A„ grow linearly proportional to the 
scaling factor x. In the limit of a large x only the diagonal 
term with m = n survives in the sum in Eq. (|19[) , and we 
have a local generator L that only retains the part of J z 
diagonal in the eigenbasis of K, 



n 

Numerically, we correspondingly see that the maximum 



precision converges to a constant, and with N — > oo ap- 
parently to a constant fraction of the Heisenberg limit. 

The remaining limits we consider are for varying e with 
the other parameters held constant. In the case e — > oo 
the other terms in the Hamiltonian become insignificant 
perturbations, and the Heisenberg limit as for K = eJ z 
is reached. In the contrary limit e — > the measurement 
precision tends to a nonzero constant. The way it works 
in the lowest order in e is seen from Eq. (p~9|) : Calcu- 
late the vectors \n) and eigenvalues A„ simply by setting 
e = in the operator K, find the local generator L, and 
obtain the maximal Fisher information from the range 
of its eigenvalues. This limit is of some interest as de- 
tection and measurements of small energy differences e is 
conceivably a frequent task. It is, however, not worth the 
space to draw a new figure, as the case e ~ is effectively 
realized whenever e < J and e < If plotted with the 
same axes as in Fig. [5J the e = results would be vir- 
tually indistinguishable from the e = 1 results already 
shown. 



V. CONCLUDING REMARKS 

It should be understood that we approach the Heisen- 
berg limit purely from the perspective of principles. Ex- 
cept for one incidental remark we pay no heed to the 
question of preparing the optimal probe state, nor to 
practical measurement strategies, nor to the problems 
with imperfect real experiments. With these caveats, 
we have studied the combined effects of nonlinearity and 
"arm-to-arm" tunneling on the best possible measure- 
ment precision of an interferometer. We use measure- 
ments of the difference of the energies of the potential 
wells in a double- well trap for bosonic atoms as a tangible 
example. We have adapted time independent perturba- 
tion theory to prove that the Heisenberg limit cannot be 
beat indirectly as a result of the nonlinearity due to atom- 
atom interactions, and to study numerically the behavior 
of the measurement precision with varying problem pa- 
rameters. Many of the limiting cases are also easy to 
understand on the basis of the perturbation theory. As 
an interesting aside, we have noted that in the limit of 
very strong attractive atom-atom interactions the input 
state required for Heisenberg limit precision is in fact a 
zero-temperature ground state. 
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